Bound States in the One-dimensional Hubbard Model 
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The Bethe Ansatz equations for the one-dimensional Hubbard model are reexamined. A new 
procedure is introduced to properly include bound states. The corrected equations lead to new 
elementary excitations away from half-filling. 
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The one-dimensional Hubbard hamiltonian, 



+ h.c. 



(1) 



describes fermions moving on a ring of length L with an 
on-site interaction U. The model was solved by Lieb and 
Wu Housing Yang's generalized Bethe Ansatz method 
(BA) ||]. They derived the following equations, describ- 
ing an eigenstate of (|^) in the sector with N particles and 
total spin S = ^{N ~ 2M), 
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A-v — sin k-i 



sinfcj + ij 
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where u — U/t. The Af "spin momenta", Xj, and the N 
" charge momenta" , kj , describe the many-body correla- 
tions of this state. Once these parameters are determined 
from the BA equations the energy and momentum of the 
state are, E — —2tJ2j coskj, P — kj. See ref [|| for 
a detailed derivation. 

We shall argue that the equations as they stand are 
incomplete. The problem arises when one considers com- 
plex momenta solutions which describe bound states and 
double occupancy. In the standard approach these 
momenta, k^ = q± are driven, in the infinite volume 
limit, L — > cxD, to "string positions" 



sin k^ — Ai ^ i 



(4) 



where e = 0{e~~^) and A; is identified with an element 
in the set of spin momenta {A^}. This expression is then 
inserted into (||,||). Simple manipulations lead to. 
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The kj, j — 1 . . . iV" denote the real momenta, kf = 
qi ± I ~ I . . . N'', are paired with the Aj and A-y, 7 = 
1 . . . Af", runs over the set of unpaired spin-parameters, 
TV = TV" + 2N'' and M = A/" + N''. The phase e''<" is 
defined as 



A; — sin k 



Ai - sin/c, + ij 



We show now that the "string Lieb-Wu" equations, 
(I'HjZil)' cannot describe all (highest weight) states in 
the infinite volume limit. In other words, the infinite 
volume limit cannot be obtained as i — > cx) limit of 
the solutions to (§,||) for finite L. The reason is, that 
equations (|^), which do not contain L explicitly, become 
redundant for L — s- 00. Maintaining them, as is com- 
monly done, cuts out a set of eigenstates which do exist 
in the infinite system. To show this, we use the algebraic 
Bethe ansatz; assume a pair of complex conjugated mo- 
menta A:"*", k~ related to the spin momentum A by (jj). An 
eigenstate of the corresponding inhomogeneous transfer 
matrix Z{p,) with (arbitrary) spectral parameter /i then 
has the form (see e.g. @): \ip) = 0^=1 B{Xj)B{A)\uj) for 
total spin S = 5 (TV" - 2Af"), where the M" creation op- 
erators B(Xry), acting on the ferromagnetic vacuum \lu), 
create A^" j-spins, and B{A) creates the j-spin of the 
bound state. By explicit calculation we find that B{A) 

acting on |a;) diverges as for L ^ 00, whereas the 

_B's not associated with the complex pair do not diverge. 
Normalizing j?/') by multiplication with e~~^ yields ac- 
cordingly an exponential suppression of vectors of the 
form Yls B{Xs)\uj) which do not contain B{A) among the 
B{Xs)- Now the equation (j^) for A is necessary to cancel 
an "unwanted term" in the eigenvalue equation for \il>): 
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(Z(/x)-£;(m))|V) = 

M" M" M" 
7 i5#7 7 

It ensures, in particular, that aQ = 0. But, as argued 
previously, the vector Y[j B{\^)B{fi)\u!) is projected out 
in the infinite volume limit (exponentially suppressed for 
finite L), and therefore (Q) is not necessary for to be 
an eigenvector of Z[fi) for i — » oo. That means that only 
a subset of all states allowed in the infinite volume can 
be generated by the L — > oo limit of solutions to (||,||) for 
finite L. 

We proceed to introduce a procedure that generates all 
(highest weight) eigenstates of the hamiltonian in the in- 
finite volume limit. We begin by deriving the elementary 
bound states of the hamiltonian, and constructing the S- 
matrices between them and the unbound particles. The 
S-matrices will be employed to derive the consistent BA 
equations. The normalizable two-particle bound state 
solutions of the hamiltonian have the form, 



gig(ni+n2)g-C(g)|ni-n2| 



(9) 



where the parameters q and ^ are related by, sinh^((7) — 
~4cos9- ^ > 0, the range for q in the interval [— tt, tt] 
is: 1^1 > f for J7 > and \q\ < | for [/ < 0. In the 
former case the real parameter q (which is half the quan- 
tum mechanical momentum: p = 2q) lies near the center 
of the Brillouin zone and in the latter near the edges. 

The bound states are necessarily spin singlets, and 
their energy is E ^ — 4i cos g cosh ^ = [/coth^(g). Thus, 
it depends directly on [/, and is negative in the attractive 
case {U < 0) and positive in the repulsive case. Moreover 
there is a gap in the spectrum: E{q) > \U\ for U > 0, 
and E{q) < -\U\ for U < 0. The hmiting values ±\U\ 
correspond to q ±{tt/2), with ^{q) tending to infinity 
at these points. 

Typically, 2-particle eigenstates of the model have the 
form (denote ni^2 = ni — n2 and A the antisymmetrizer), 

F{ni,n2) = 

= ^e*fci«ie^fc2"^ [Aa,a,0{n2s) + Si2Aa,a,0{ni^2)]. (10) 

with 5i2 = Si2 being the S-matrix between unboundpar- 
ticle 1 and 2 carrying real momenta ki and fc2- It reads. 



^12 



(fci, h) 



sin fci 



sinfe + i|-Pi2 



sin ki — sin A:2 -f i-^ 



(11) 



We can view the bound state as corresponding to com- 
plex momenta fci = fc~ = g — i^, fc2 = fc"*" = g + i^, with 

/ 2 \l/2 

sinA:=t = 0(g) T if, where (^(q) = sing \\ -f jg^^j . 

Thus 5*12 vanishes for the (singlet) bound state, while 6*21 
has a pole. Eq ( p^ ) therefore gives correctly the wave- 
function (^ after antisymmetrization because the "for- 
bidden" function e^" in interval n > is projected out. 



One sees that a subset of S-matrices become undefined, in 
this case S'21, within the Bethe ansatz for single-particle 
states for L — > 00 ||] . To avoid these singular S-matrices 
one has to treat the bound pair as a composite with a 
well defined S-matrix S"^^ with unbound states as well as 
with other bound states, S^^ . To compute 5*"^ we write. 
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ub 
(23) 



{k, k'^)Si2 {k, k ) 
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and note that the spin space of the three particles, 
V, is restricted: 2 and 3 form a singlet in the region 
[123] denoting the ordering {ni < 712 < 713). It is 
V — Vi ^ V^g*"^''^* J constituting a two dimensional sub- 
space of the eight dimensional spin space of the three 
particles. In this space -P23^[i23] = ^^[123], so that 
P13P12 = ^'i2-P23 = -Pi2- Further P12 = 1 - P13 when 
acting on V. In general these operators do not leave V 
invariant, and the singlet state of particles 2 and 3 will 
be destroyed upon scattering with particle 1 . It is due to 
a non trivial cancellation of terms for momenta fc+, /c~ 
satisfying the pole condition for 5*23*, that £'"('233 indeed 
leaves V invariant, acting as a phase on the wavefunction. 
We find, 



'^1(23) 1) 



In an analogous way, 



sin k — (j){q) — i\ 
sin k — (j){q) + i\ 



Hqa) - (j>{qp) + i^' 



(13) 



(14) 



Proceeding in the same manner we find more elaborate 
bound states: bound states of bound states. We can 
infer the existence of these higher composites by look- 
ing for zeros, respectively, poles, in the S-matrix for 
two two-particle bound states (|lj). The zero of (|^ is 
at (t){qa) — 4>{<li3) — if, corresponding to complex g's. 



Choosing 



± ij, we find the four momenta 



of this "double" bound state, the so-called quartet 



sinfc[^2 — 4>X' , sinfcg^^ ~ 
this state with an unbound particle is, 



(2) 



hf' ±i\. The S-matrix of 



~iu{hh') 



ifc 
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In general there are bound complexes of 2m electrons, 
(m-complexes) , which correspond to a pole in the S- 
matrix of an (m — l)-complex with a simple bound pair, (a 
1-complcx). The m-complex is parameterized by an m- 
string of the form: 



+ (m-hl-2j)z- j = l. 



(16) 



The S-matrix of an unbound particle with a m-complex 
is 
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u{m) 



sin k 



sin k 



a(™) 



(17) 



and the S-matrix of an m-complex with an n-complex: 



S 



("')(") 



- ^/-i"^ + |n - - + (n + m)z 



locahzation length of all bound states tends to zero, see 
(^), periodicity in the center-of-mass coordinate of the 
composite becomes equivalent to strict periodicity. From 
continuity follows then that the total number of states is 
4^ for any finite u. 

J Diagonalizing the transfer matrices leads to the follow- 
ing set of equations, 
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min(m,n) — 1 / Am) 



n 



{\n-m\ + %)i\ 



Using these S-matrices we construct, in the usual man- 
ner (see e.g. [||), eigenstates of the hamiltonian charac- 
terized by: A'^", the number of unbound particles, and 
7V^, m = 1 . . . cxD, the number of m-complexes. These are 
indeed eigenstates since they derive from S-matrices sat- 
isfying the Yang-Baxter relation. The antisymmetriza- 
tion of the wavefunction is trivial within a complex, and 
amounts to a change of labels between the constituents 
of the complexes. 

We now impose periodic boundary conditions and de- 
rive the BA equations. Consider a system with N = 
iV"-|-2iV^ electrons, where TV" particles are in plane wave 
states given by momenta kj, j = 1 . . . iV" and 2N^ parti- 
cles are in bound states, N'' = mN^. The total spin 
S is given by S" = i(7V" - 2M"), where M" denotes the 
number of down spins among the iV" unbound electrons. 
Imposing periodic boundary conditions on the unbound 
particle j we are led to the diagonalization of the matrix 
Zj which takes it around the ring. 



QUU 



QUU 



ii(m) 



.S 



ii(m) 



(18) 



with eigenvalue Zj = e^^^^. Similarly, taking an m- 
complex around the ring leads to the matrix (actually 
a phase) Z^a,7n} 



ti(m) 



{b,n] 



(ri)(m) 



(19) 



where the index (6, n) runs over the set of all complexes 
present (omitting (a,m)). The corresponding eigenvalue 
is Z(^a,m) = exp (iL[— 25ftarcsin(0a -I- mi^) -I- (m -I- l)7r]). 
The modified boundary condition for the bound com- 
plexes allows for a total number of states, 4^. This can 
be seen as follows: The number of states in the system 
does not depend on u, but since in the limit u ^ oo the 
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Xs — sin kj 



Xs — sin kj 



(m)t 



-| — r A-y \^ Z"2" T — r A--y SlTl kj Z 

|i X^-Xs + q " \}^ A^-sinfc, +if 



u{ra) 
kj 0, 



(20) 
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(22) 



The equations ( [20[|2l| , |22| ), to which we refer to as Bethe 
Ansatz equations for Composites (BAC), determine the 
eigenstates of the Hubbard hamiltonian in the presence 
of bound states. They are weaker than the "String Lieb- 
Wu equations" . The latter include additional constraints 
which lead to over-determination; a quick way of estab- 
lishing it is to carry out a counting of the states described 
by the BAC equations. This can be done along the lines 
of 1^,^, and one finds exactly 4^ states, proving a forte- 
riori that the stronger set (|[^ is over-determined.^ 

The BAC equations modify the structure of the up- 
per Hubbard band. We shall show now that in addition 
to the holon-antiholon excitation which have linear dis- 
persion at low energies a new excitation appears with a 
nonlinear dispersion. Consider a single bound pair above 
the ground state of the repulsive Hubbard model. To 
keep the number of particles N unchanged two holes are 
placed in the sea of unbound particles. Thus, A^" = N — 2 
and M" — N/2 — 1 in the notation above. Equations 
(|^,|l|,|2|) become, 

Lkj = 9i{smkj — A^) + di{smkj — (f>{q)) + 2iTnj 
s 

^ 0i(sinfcj - A-^) = ^92{Xs - X^) + 2ttI^ 



(5=1 



JV-2 



2qL = ^ 9i{(j}{q) - sinfcj) -f 27rJ 
i=i 



^Takahashi @ obtained the set (^||) from the set (|5|-|) 
neglecting the terms of order e and eliminating the (^;'s by 
plugging (7) into (8). He assumed the same number of solu- 
tions as follows from BAC. In this way he was able to derive 
the correct thermodynamics. 
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where 9n{x) = —2 tan (f :^)- The range for the quan- 
tum numbers nj is -N/2 < rij < N/2 - 1 for 
even, and -{N - l)/2 < < {N - l)/2 for M" 
odd. The range between -((iV - 2) - - l)/2 
and +{{N - 2) - M" - l)/2. The nj sequence con- 
tains two unfilled slots, holes, as the actual number of 
unbound fc-momenta is N — 2. We denote the hole po- 
sitions by ki and k!^. The Xy sequence does not contain 
holes in the absence of spin excitations. J, the quantum 
number associated with the bound state, is an integer if 
L — N is even, and a half-odd integer if it is odd. We 
assume L, N even in the following. To find the limiting 
values for J, we consider the boundaries of the allowed 
range for q: 7r/2 < <? < tt, (resp. — tt < q < — 7r/2). 
We may treat the range for q as simply connected by 
shifting -tt/2 to 37r/2. We find, L-f>J>fso 
that the associated mode becomes a full fledged exci- 
tation, independent of the two holons. At half-filling: 
{N/2) < J < {N/2), which leads to maximal restric- 
tion of the phase space for the bound state. To com- 
pute the spectrum we proceed in the usual manner [0J|] 
converting the BA equations to integral equations, for 
notation see One needs to calculate p{k) the den- 
sity of the kj momenta, which consists of several terms, 
pik) = po{k) + L-l[pi(fc; k1) + pi(fc; fct) + p\k; q)]. 

The first term describes the kj distribution in the 
ground state, pi{k;kj), j = 1,2 are the holon contribu- 
tions, and p^{k;q) is the bound-state contribution. We 
write here only the integral equation for p'^{k; q), 

p^{k)-cosk- r dk'R f-(sinfc-sinfc')^ P^(fc') 
Uj-Q \u J 

= cosfc Ki{smk — (t){q)). 

Here Q is the k Fermi-level determined by the require- 
ment p{k) dk = N''/L. 

We find that the total excitation energy and momen- 
tum consist of contributions from the two holons and 
from the boimd state, 

AE^E-Eo = e'^kf^) + e'\k'i) + e\q) 

AP = P-Po^ p''{k'^)+p\k!^)+p''{q) (23) 



holon parameters fcj*, k!^. However, away from half- filling 
the bound state becomes an independent excitation with 
momentum, p''{q) = 2q — J^Qdkpo{k)9i{(j){q) — sinfc). 
As q runs over the allowed range 7r/2 < g < tt (resp. 
— TT < q < —tt/2) , p^ varies in the range (tt -|- A;F,27r) 
(resp. (— 27r, — tt — kp). This corresponds to a symmet- 
ric band around zero between —\tt — kpl and n — kp. 
It follows that at half-filling p^ = 0, the phase volume 
vanishes. Starting directly from this case and setting 
J — 0, q becomes fixed in terms of the hole parameters 
fcj' and k^: (j){q) = ^(sinfcj' + sink!^) [|. This state- 
ment, however, is empty as the physical parameter of the 
excitation, the energy and momentum, are independent 
of the unphysical parameter q: = U and p^ = 0. 
Away from half filling the dispersion of the new mode is 
given parametrically by f^{q) and p^{q). The solution of 
the equation for p** has to be performed numerically in 
general. However, at the band edges, p^ = ±(7r — kp), 
analytic results are accessible. In particular, the velocity, 

= (de^/dp'')|pb^±(jr-fcF) = 0. It follows that this ex- 
citation is quite different from the gapless particle-hole 
excitations, as the nonzero curvature of the dispersion 
at the band edges sets an intrinsic scale. Because these 
excitations do not appear in the low energy regime they 
do not modify the transition to the known continuum 
limit of the model (the Luttinger liquid but the 

parameters of the model need to be further studied, [ p^ . 

Similar consideration lead to a modification of the 
(gapped) spin excitation sector in the attractive case. 
Again, a new dispersive spin singlet mode makes its ap- 
pearance in addition to the spinous |^ . Also other models 
where complex momenta appear, such as the multichan- 
nel Kondo, the Anderson or the t — J model, need to be 
further examined |p^ . 
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with direct, and a backflow contribution. The holon exci- 
tation is well known |^,||. In particular, e^{k'^) vanishes 
linearly as A:'* — > Q. This corresponds to the holon mo- 
mentum, Pj— /q ^ dkpo{k), tending to the charge Fermi- 
momentum ±Kp = ±2fci? = ±7r(7V/L). The bound state 
energy is computed from. 



e'iq) = £{q) - 2t 



Q rQ 

Akp^{k;q)cosk — p I dkp''{k;q) 



with £{q) — — 4t cos(7Cosh(^((7)), and p — dE/dN, being 
the chemical potential 0,^. At half-filling e''{q) = U and 
the dispersion of the total excitation depends only on the 
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